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Abstract 

We consider scalar field theories in dimensions lower than four in 
the context of the Wegner-Houghton renormalization group equations 
(WHRG). The renormalized trajectory makes a non-perturbative in- 
terpolation between the ultraviolet and the infrared scaling regimes. 
Strong indication is found that in two dimensions and below the mod- 
els with polynomial interaction are always non-perturbative in the 
infrared scaling regime. Finally we check that these results do not 
depend on the regularization and we develop a lattice version of the 
WHRG in two dimensions. 
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1 Introduction 



Asymptotically free models gained importance both in High Energy and in 
Condensed Matter Physics. This is obvious for the latter, where the low 
dimensional phenomena invoke effective theories below the upper critical 
dimension, d — 4, which are super-renormalizable and their Hamiltonians 
contain asymptotically free coupling constants only. In the former case the 
the asymptotically free models can provide structureless high energy physics 
where the cutoff can be pushed away at will. What kind of non-perturbative 
mechanisms do we encounter as we follow the renormalization group flow of 
an asymptotically free model? To distinguish different finite energy mech- 
anisms we recall that unless the model possesses unbroken scale invariance 
there are always two scaling regimes, an ultraviolet and an infrared one sep- 
arated by a crossover at the characteristic scale of the model. If there is a di- 
mensional coupling constant, say a mass parameter m, in the Lagrangian then 
the UV and IR scaling regimes are separated by a crossover at p = p cr ~ m 
because m 2 /p 2 or p 2 /m 2 is treated as a small quantity in the UV or the IR 
side, respectively. The IR scaling laws are trivial because the radiative cor- 
rections to the evolution are suppressed by p 2 /m 2 and what is left is the scale 
dependence governed by the canonical dimensions of the coupling constants. 
If there are no dimensional coupling constants in the Lagrangian, i.e. the 
model possesses classical scale invariance then a crossover can be identified 
where one of the coupling constants reach the value 1, g{p C r) = 1- The ex- 
pression for the dimensionless running coupling constant g(p) must contain 
a dimensional parameter, usually called the A-parameter and p cr ~ A^J] 

Non-perturbative effects may originate in the IR scaling regimes. The 
classically scale invariant models at the lower critical dimension, such as the 
two-dimensional sigma and the Gross-Neveu model do not support long range 
order. Thus the IR scaling laws must contain a dynamically generated mass 
scale, an effect generated by the infrared or collinear divergences. 

The asymptotically free coupling constants may generate non-perturbative 
effects in the UV scaling regime, as well, due to their growth. The question is 

*It is worthwhile noting that one can introduce the A-parameter even in the presence of 
any parameter with positive mass dimension in the Lagrangian so long as the IR dynamics 
remains stable when this latter is removed. In fact, let us write this parameter as m where 
m defined in this manner is the mass scale of the model. The running coupling constant 
depends on the ratio m/p, g(p) — f(m/p) where the limit m — > is convergent. In other 
words, the evolution of the classically scale invariant system is recovered as p — > oo. 
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whether the scale parameter A^t or the mass m is reached first as we lower the 
cutoff. The system remains perturbative for m > because the IR scaling 
laws cut off the growth of the asymptotically free coupling constants before 
they would reach a dangerously large value. The dynamics of the modes 
p < A^4 becomes non-perturbative when A^4 > m. The four-dimensional 
Yang-Mills models develop linearly rising potential between static external 
charges. This is believed to happen due to the large value of the asymptoti- 
cally free coupling constant at the ultraviolet side of the crossover, indicated 
by the infrared Landau pole in perturbative QCD. The quark masses are 
supposed to play no important role in the vacuum structure and the chiral 
limit m — > is assumed to be safe and convergent, though non-perturbative. 

The infrared singularities are easier to isolate than the non-perturbative 
effects arising at the IR end of the UV scaling regime. It has already been 
achieved for non-asymptotically free models |]J and superrenormalizable the- 
ories ||. The more complex non-perturbative features should arise at the 
ultraviolet side of the crossover where no asymptotic analysis is available. 

The goal of the present paper is the study of the scaling laws and their 
consequences at the IR end of the UV scaling regime. According to our best 
knowledge this source of the non-perturbative effects in asymptotically free 
theories has not been studied before in a systematic manner. The problem 
is rather involved and demanding, we restrict ourselves to outline only the 
way the more detailed analysis should be done. In order to defuse the in- 
frared problem we turn to the scalar model with polynomial interaction in 
the phase where no spontaneous symmetry breaking occurs. The term <p n 
in the Lagrangian of the d- dimensional scalar model is asymptotically free 
for n < 2d/(d — 2), so we have to trace the evolution of a large number of 
coupling constants at low dimensions. We present numerical evidences that 
the polinomial interactions at the lower critical dimension or below where in- 
finitely many operators are relevant in the UV scaling regime always become 
strong, non-perturbative. 

The study of the asymptotically free three-dimensional 6 model can bring 
some light into the IR Landau pole problem of four- dimensional gauge the- 
ories by considering the role of the non-renormalizable couplings in the IR 
extrapolation of the UV scaling regime. We will show that they affect consid- 
erably the evolution of the relevant couplings. This points to the possibility 
that the non-renormalizable couplings could suppress the Landau pole. The 
full solution of the theory where all possible coupling constants are followed 
by the renormalization group equation should necessarily yield a non-singular 
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renormalization group flow so long as the theory possesses local interactions. 
This suggests that the introduction of the hadronic composite operators in 
QCD could defuse its Landau pole problem [[J. 

Another issue addressed in this paper is the manner in which the influence 
of the non-renormalizable operators on the dynamics is suppressed and the 
universal physics is reached as the cutoff is lowered. The usual argument 
based on the linearization of the blocking relation is not obviously applicable 
due to the presence of infinitely many relevant operators. This question is 
discussed in the framework of lattice regulation by tailoring the Wegner- 
Houghton equation for lattice regulated models. As a result, the approach 
to the low energy physics can be compared for the momentum space cutoff 
regularization in the continuum and for the lattice regularization, and the 
universality, the regulator independence, can be established well beyond the 
linearized approximation of the blocking relations. 

The organization of this paper is the following. The infinitesimal renor- 
malization group step, the Wegner-Houghton equation, is described in Sec- 
tion II. The running coupling constants are introduced and their evolution 
equations are given in Section III for scalar models. The asymptotic UV evo- 
lution is discussed in Section IV for d = 2, 3 and 4. Section V is devoted to 
the demonstration of the difficulties in finding a perturbative model in d = 2. 
The Wegner-Houghton equation is derived in the lattice regularization and 
its solution is presented in Section VI. Section VII is for the conclusions. 

2 The Wegner-Houghton equation 

There are different ways the mixing of a large number of operators can be 
traced down. The Wegner-Houghton equation [Q, which we use in the local 
potential approximation in this work, is the simplest implementation of the 
Kadanoff- Wilson blocking || in the momentum space and produces the cut- 
off dependence of the bare action along the renormalized trajectory. Other 
methods work with the effective action where the infrared cutoff dependence 
is sought |J. Different schemes should agree in the infrared limit where few 
long wavelength modes are left only in the system. We shall make two ap- 
proximations in computing the blocked action, the truncation of the gradient 
expansion at the leading order, the local potential approximation ]7[] , and the 
truncation of the Taylor expansion of the local potential in the field variable 
|| . The higher order terms of the gradient expansion are non-renormalizable 
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according to the power counting. We believe that these coupling constants 
which are irrelevant in the ultraviolet scaling regime do not modify our qual- 
itative conclusion. 

The Wegner-Houghton (WH) equation [|J describes the evolution of the 
effective action as the cutoff is lowered. As mentioned in the Introduction 
we shall consider a scalar model with an intrinsic mass scale which allows 
to clearly distinguish the UV and IR scaling regimes. We derive the WH 
equation for a scalar field theory by using a sharp momentum space cutoff |3|, 
[|: we will call this regularization procedure the continuum regularization. 
In Section |6| we consider an alternative, the lattice regularization. 

Denote the bare action by [4>] , where k is the UV cutoff. Then, accord- 
ing to the usual Wilson-Kadanoff procedure, 

e -\s k ,m = f V(f) > e -is k \^'] (i) 



where k' < k in the Euclidean space-time. The Fourier transform of the fields 
4>(x) and (p'(x) is non- vanishing only for p < k' and k! < p < k, respectively. 
The right hand side is evaluated by means of the loop expansion, so that 
Eq. (P gives 

SM = S k {0 + <f>' ] + 5tr log 5 2 S + 0(h 2 ), (2) 



where 



fiSlXiy)= ™&+M, (3) 



S4/(x)5(j/(y) 

and the saddle point, <p' , is defined by the extremum condition 



ss k [<i> + <f>r ] 



0, (4) 



in which the infrared background field, <f)(x), is held fixed. It can be proved 
that the saddle point is trivial, 4>' Q = 0, as long as the matrix 5 2 S(x, y) is 
invertible and the IR background field is homogeneous, <p(x) = $. 

Now, each successive loop integral in the n-loop contributions which are 
not explicitly written in Eq. brings a suppression factor 

kd ~ k,i ( k -±) (5) 



k' d V k 
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due to the integration volume in the momentum space. Thus Sk/k = (k — 
k')/k appears as a new small parameter which suppresses the higher loop 
contributions in the blocking relation and the "exact" functional differential 
equation obtained in the limit 5k — > includes the one-loop contribution 
only. But we should bear in mind that the loop expansion had to be used 
at the initial stage of the derivation so the resulting "exact" equation might 
be unreliable in the strong coupling situation. All we know is that the loop 
corrections to the evolution equation obtained in the one loop level are van- 
ishing. 

We will use the gradient expansion for the action, 

CO „ 

S[<f>]=J2 / d rf xf/ n (0(x),9 2n ), (6) 

n=0 ^ 

where U n is an homogeneous function of order 2n in the derivative. In the 
leading order of this expansion, the so-called local potential approximation, 
we have 



S\6]= / d d x 



(7) 



and furthermore the simplification Z(<f>) = 1 will be used to derive a simple 
differential equation for the potential U. This local potential will be then 
the only function characterizing the action. If we use now a homogeneous 
infrared background field, <p(x) = $, we obtain from Eqs. (Q) and (0) an 
equation for the local potential Uk(<f> = 3>): 

U k -6k($) = U k ($) + itrlogp + + 0(6k 2 ), (8) 

where we have introduced the notation 

and the trace is taken in the subspace of the eliminated modes. We can 
explicitly write the trace in momentum space and get 

tW$) = u k @) + lj-0r d Mp 2 + Kin + 0{8#), (io) 
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where the integration extends over the shell k — 5k < p < k. In the limit 
5k — > one then finds the differential equation 

^ ($) = -^F log[fc2 + ^ ($)] ' 

where fl^ denotes the (i-dimensional solid angle 

n " = m (12) 

The Wegner-Houghton equation (|ll|) represents the one-loop resummed 
mixing of the coupling constants of the potential Uk = ^2 n (g n /n\)$> n . In 
fact, an expansion of the logarithm in the second derivative of the potential 
gives 

% w = -w^U + wJ ' (13) 

up to a field independent constant. This is the usual one loop resummation of 
the effective potential |l(J except that the loop momentum is now restricted 



to the subspace of the modes to be eliminated. Actually, the fact that the 
r.h.s. includes the running potential Uk{<&) rather than the bare one, £/a($), 
indicates that the contributions of the successive eliminations of the degrees 
of freedom are piled up during the integration of the differential equation and 
the solution of the renormalization group equation resums the perturbation 
series. The solution of the differential equation interpolates between the bare 
and the effective potential as k is lowered from the original cutoff A to zero. 

Finally, let us note that the derivation of Eq. (|TTD shows that the restoring 
force for the fluctuations into the equilibrium is proportional to the argument 
of the logarithm function. Thus a nontrivial saddle point should be used when 

k 2 + < 0. (14) 

3 Evolution of the coupling constants 

Our effective action is 



S, = / d d x 



(15) 
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and the initial condition for the evolution equation is given at k = A. The 
potential £4(<&) is assumed to be polynomial so it is expanded as 



N 1 

=V^?n($ )($-$0) n . (16) 



n=0 



We study the model in the symmetric phase, where the saddle point is trivial, 
$o = 0. The polynomial structure of the potential is consistent because we 
avoid the singularity of Eq. (fLlf) at k 2 cx = —Ujl a (<&) (recall Eq. (|14])) which 
occurs in a region around $ = in the symmetry broken phase ||. 

Taking the n-th derivative of Eq. fll6|) at $ = we obtain the coupling 
constant g n {k), and define the corresponding beta function by 

9n(k) = ^£4($)|*=o, 

«• = % 9 " W = (17) 



By taking the successive derivatives of Eq. fllTD , we obtain 

2(2tt) 

where 



^ n = ~ /r> \ ri 'Pn(G2, • • • , G n +2), (18) 



and 



C„ = (19) 
k + 92 

K = ^log[P + ^($)] (20) 
is a polynom of order n in the variables Gj, j — 2, . . . , n + 2, 

Pi = G 3 , 

V2 = G4 — G\, 



V3 — G5 — 3G3G4 + 2G3, 

p 4 = G 6 - 4G5G3 - 3G 2 4 + 12G 2 3 G 4 - 6Gj, (21) 
V 5 = G 7 - 5G 6 G 3 - IOG5G4 + 20G 5 Gl + 30G^G 3 - Q0G 4 G 3 3 + 2AG 5 3 , 
V 6 = G 8 - QG 7 G 3 - 15G 6 G 4 - 10G 2 5 + 30G 6 G 2 3 + 120G 5 G 4 G 3 + 30G 3 4 
-120G 5 G 3 - 270GlG 2 3 + 360G 4 G% - 120G|. 
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The coupling constants defined through Eq. flTUD are dimensional param- 
eters (the field variable has dimension (d — 2)/2). However, the correspond- 
ing dimensionless parameters have more physical sense. We obtain them in 
the following way: 

g n (k) = k~< l ~^ d g n {k) = (Ak)- n ^- d g n (k), (22) 

where now k runs from 1 to 0. Their beta functions are 

~g n + k~ n ^- d (3 n , (23) 

where the first and the second term stands for the tree-level and the loop 
corrections, respectively. One can see that the super-renormalizable coupling 
constants follow asymptotically free scaling law at the tree level. 
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4 UV scaling laws and their extensions 

One can distinguish an ultraviolet and an infrared scaling regime, for k 2 ^> 
\m 2 (k)\ and for k 2 <C \m 2 (k)\, respectively. In the UV regime the scale 
dependence comes dominantly from the k 2 term of the propagator, see the 
denominator of Eq. ([H]); the /c-dependence is generated by the phase factor 
k d in the IR regime where k 2 could be neglected in the inverse propagator. 
We will begin at the UV scale with the usual 4 potential (#2 = rn 2 ) 

V A (<P) = \m 2 4> 2 + 1^ 4 4 + ^ 6 6 , (24) 

and see how the different couplings are generated when we move towards the 
IR regime. 

One ignores the #2 term in the denominator of Eq. ([19]) in the asymptotic 
UV regime and finds 



d#2 


fid 


dA; 


2(27r) d 




Qd k 


dk 


2(2vr) d 


dg 6 


n d 


dk 


2(27v) d 



k a - 6 gi, 



d-3 I 3 o 



-j^9i-g&U (25) 



3g*k d - 5 ( ^U 2 - 5g< 
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where in the last equation we omitted the contribution of g s . Consider the 
usual strategy in which the coupling constants g n are neglected for n > 4 
and the resulting equation is easy to integrate, 

j_ = J _ an, (i - (I) 4 " 1 ) 

94(A) 94 (A) + (4-,i)2(2,r)'<A 4 -<'' (26) 

This expression agrees with the result of the minimal subtraction (MS), a 
scheme which proved to be specially convenient in the ultraviolet scaling 
regime. It is based on the analytical continuation of the loop integrals in 
the ultraviolet domain so the resulting beta functions are mass independent, 
i.e. the terms 0(g2/k 2 ) are neglected. When extrapolating to the infrared 
regime we find erroneously the mass independent result ~ k A ~ d (g 4 ~ 
log A; in d — 4), g± tends to zero as k — > 0. This can be understood by 
inspecting ( P3|) where we find large positive values in the infrared for g 4 7^ 
(this conclusion remains valid for finite g§, as well). When the mass term is 
retained the beta function assumes the correct behavior and becomes 0(k d ) 
in the infrared. Note that the term with g§ acts in the opposite manner than 
g 4 , c.f. the different signs in the right hand side of fl2"5|), so that it can change 
the evolution considerably. 

It is instructive to look into the evolution of the dimensionless coupling 
constant, 

j_ (r y"-('-<r) 

Hk) 94(A) ~" (4 - d)2(2ir)<i ■ 1 ' 

For 4 — d > the one-loop uo = 4 — d universal critical exponent is reached 
for k values sufficiently below the cutoff where k/A w 0. The latter condition 
is needed to get rid of the non-universal cutoff effects. The scaling changes 
qualitatively as d — > 4 because the non-universal k w A behavior is spread 
over the whole k range due to the smallness of 4 — d. This is what happens 
in the expansion 

l-(£f"-(rf-4)h4, (28) 

employed in the dimensional regularization scheme. This generalizes to any 
dimension: the marginal coupling constant follows the scaling law extended 
from the non-universal cutoff regime. 
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The evolution for g 2 is of the form 

94k d ~ 3 . (29) 



dg 2 _ tt d , 3 



dk 2(2vr) 
It predicts 



2 



dk 



k, g% ~ k 2 + const (30) 



in lack of any dimensional constant. 

In the IR scaling regime we neglect the k 2 term in the denominator of 
Eq. (|i9|) and using Eqs. (|18D and ([21]) , we get for g 4 (assuming again g 6 — 0), 



dc/4 _ 3fl d d _! ^ 

d* - 2(2tt)- A gl [dl) 

This evolution is much slower comparing with Eq. (^). In fact, for d > 1 we 
have now a suppression factor A; d_1 which makes the coupling to stabilize at 
the attractive IR fixed point. In the same way, we obtain for g 2 

dg 2 _ _ rd-i 94 mm 

dk ~ 2(2vr)^ g 2 ' ^ 

with a variation which is slower than that predicted by the UV scaling and 
a suppression factor for k — ► 0. 

We see that the extraction of the scaling in the limit k —>■ from the 
UV scaling laws, which is the commonly accepted practice in perturbation 
theory is incorrect when the MS scheme is used. One has to come back to 
the complete scaling laws in order to describe correctly the IR scaling. 

Four-dimensional asymptotically free gauge models present an infrared 
Landau pole in perturbation theory. But this behavior results from the ex- 
trapolation from the UV scaling laws. As we have just remarked, the IR limit 
of the UV regime is not correct in general, and the mass term can change 
considerably the actual behavior in the IR. Moreover, at the IR side of the 
UV regime, there are nonlinear effects that make important the contribution 
of the irrelevant (non-renormalizable) couplings (see Fig. 3 (b), commented 
in the next Section), therefore even the IR limit of the UV scaling can be 
influenced by these couplings. These ideas have been considered qualitatively 
in the previous paragraphs, after inspection of Eq. (p5|). Let us now examine 
them more quantitatively. 
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We take as an example of asymptotically free model the scalar theory in 
three dimensions. We know from the epsilon-expansion result [llj that below 
four dimensions, the A0 4 theory does not present an infrared Landau pole, 
but a fixed point located at 

~ 167T 2 . , 

A* = — e, (33) 

at order e = 4 — d. This is the Wilson-Fisher fixed point. For finite e, for 
example in two and three dimensions, we can get this one-loop result from 
our beta functions of the dimensionless coupling constants obtained from 
Eqs. ( p3|) and (JfJ) in the asymptotic UV regime (that is, ignoring g 2 ), which 
for the #4 coupling give 

P, = -(A-d)~g 4 + 3~gl^y d , (34) 



giving the IR fixed point 



„, (4 - d)2(2K) d 

9 * = sn d • (35) 

Restricting ourselves to the d = 3 case, we find 

4tt 2 

~g\ = — ~ 13.15947. (36) 

However, in three dimensions g§ is a marginal coupling, and it can be gen- 
erated by the RG flow, modifying the position of the fixed point (|36|) . If we 
include c/q in our analysis, we get the beta functions 

P* = ^4(10^ 4 2 -5^ 6 ). (37) 



It is immediate to see that the zeros of these beta functions are at the point 

g* 4 = 4tt 2 ~ 39.4784, 

~gl = 32vr 4 ~ 3117.091. (38) 
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But now let us consider the inclusion of a non-renormalizable coupling, g s . 
The beta functions for the (g 4 , g G , g 8 ) model are 

^ = + 4^ ~ 4^ 6 ' 

— 3 _ 1 
h = -^94 (Wgl - 5g 6 ) - —g 8 , (39) 



Ps = 98 ~ ~r~~2 (~ 9 °9t + 60gjg 6 - 5g% - Ag 4 g 8 ) 



g* = ( 195 + V29625) ~ 69.01 



7_ 

4V 2 

From them, one obtains the fixed point 

4vr 2 
210 

~gl = -4n 2 ~gl + 3~gf ^ 11563, (40) 
g* s = 15gf - 60tt^ 2 ~ 2.11 • 10 6 . 

These values are also obtained in the numerical integration of Eqs. (|39|), 
independently of the initial values for the different couplings (if they are 
different from zero, which corresponds to the Gaussian fixed point). 

To assess the importance of this result, the difference between the physics 
around the fixed points ( |38D and (|40|) , recall that the modification of the ir- 
relevant operator set at the cutoff influences the overall scale of the model. 
Thus one has to consider dimensionless quantities in comparing the two cou- 
pling constant regions. The most obvious candidate, the dimensionless ratio 
between the mass and the four point vertex, g2/gl, is trivially vanishing in our 
approximation. But g§ is dimensionless and its variation at the fixed points 
indicates that no adjustment of the overall scale could bring the physics of 
these two fixed points together. 

As we have expected, the non-renormalizable coupling g 8 modifies the 
position of the fixed point without changing the blocking procedure, turning 
to a situation of strong coupling dynamics in the IR. This is because the 
linearity which one assumes to ignore the irrelevant coupling constants is no 
longer valid in the strong coupling regime. As g± and g§ approach their large 
IR fixed point the linearization fails and new scaling laws are found which 
in turn generate new relevant operators ||, overlapping with </> 8 . Thus the 
strong coupling dynamics may induce a new (and artificial) IR scaling regime 
even if the UV scaling laws are extrapolated down to low energies. 

Nothing unusual happens for infinitesimal e when one stays in the vicinity 
of the Gaussian fixed point. In that case, linearity applies all along the RG 
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flow from the Gaussian fixed point to the Wilson-Fisher fixed point given by 
Eq. (0). However, as we have seen, its location is changed for finite e by 
nonlinear effects produced in the flow from one fixed point to the other. In 
the vicinity of this infrared fixed point (which, we stress again, is artificial 
in the sense that it neglects the influence of the mass term) we of course 
have again the classification of relevant and irrelevant terms, which, in the 
case of three and two dimensions, is different from the one obtained from 
power counting. For example, this IR fixed point has in two dimensions just 
the mass and the fourth order coupling as relevant parameters, while higher 
order couplings are irrelevant. 

It is well known that the poles of the fixed point action at complex values 
of the field variable make the Taylor expansion in the field unreliable || . We 
do not see any reason to reject a blocked action only because the potential 
is diverging beyond a given field strength. This kind of internal space sin- 
gularity might only indicate a maximal particle density in the system. We 
should stay only sufficiently far from this limiting value of the field variable 
when the evolution equation is truncated. We interpret the difference of the 
two fixed points as an indication of the breakdown of the simple universality 
which is based on the linearized flow equation around the UV fixed point. 

So far we considered the extrapolation of the UV scaling laws to the 
IR regime. Does the conclusion concerning the importance of the non- 
renormalizable coupling constant g$ remains valid when the true evolution 
equation, with g 2 7^ 0, is considered? The mass slows down the evolution 
but this may happen "too late" and the strong coupling effects can be found 
on the true renormalization group trajectory for small enough renormal- 
ized mass, close enough to the critical point. When the mass is large then 
crossover freezes the evolution of the coupling constants "earlier" and the 
linearization remains valid. To demonstrate this case recall that the IR limit 
of the UV regime means a fixed point for the dimensionless couplings. For 
a super-renormalizable coupling constant such as g±, which has positive di- 
mension, this would mean that the dimensional coupling goes to zero when 
k — > 0. However, we know that for a relevant coupling, the dimensionless 
quantity diverges when k — > 0, so that the dimensional coupling will take 
a finite value at the IR. This reasoning can be explicitly checked in Fig. 1, 
in which we consider the d = 3 scalar theory with just one coupling, g±. 
The white points follow the evolution of the UV regime and its extrapola- 
tion to k — 0. We observe that indeed the dimensionless coupling reaches 
the fixed point given by Eq. (|36D , while the dimensional coupling goes to 
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zero. However, if one considers the complete beta function, i.e. retaining 
#2 (black points), the behavior is the same in the UV regime, but then the 
true trajectory separates from the IR limit of this regime and enters into the 
actual IR regime, which implies a divergent dimensionless coupling at k — 0, 
and a certain finite value of the dimensional coupling, as explained above. 
One can however see numerically that this finite value is stable and almost 
does not change when one introduces more and more non-renormalizable 
couplings in the RG evolution. This is what one expects when the crossover 
captures the coupling constants and slows down their evolution in the regime 
of linearizability where the non-renormalizable couplings are unimportant. 

In the same way, it might well be that the IR Landau pole observed in 
four- dimensional gauge theories is just an artifact of a wrong IR limit or the 
truncation of the renormalized action. First, the nature of the singularity 
can change when one adds non-renormalizable couplings, and then the IR 
Landau pole would be the reflection of the insufficient functional form of the 
blocked action, and second, the true IR trajectory can be quite different from 
the IR limit of the UV scaling. 

5 Asymptotic freedom and the perturbation 
expansion 

We examine in this Section the scaling laws in dimensions d — 2,3 and 4 
from the point of view of the applicability of the perturbation expansion. 

As we have seen in the previous Section, in the RG evolution of our model 
there are two asymptotic scaling regimes, k —>■ oo and k — > 0. The latter one 
is trivial as mentioned above, because the beta functions ([18]) are suppressed 
by the factor k d and the evolution of the dimensional coupling constants 
slows down as k — > 0. The asymptotic UV scaling is however more involved. 
The super-renormalizable (relevant) and the renormalizable (marginal) cou- 
pling constants, g n with n < 2d/(d — 2) and n = 2d/(d — 2) according to the 
power counting, respectively, follow their autonomous evolution, the univer- 
sal renormalized trajectoryf]. The non-renormalizable (irrelevant) coupling 
constants "forget quickly" their initial value and take values which are gen- 
erated by the universal flow. This general trend is demonstrated by the 

2 Ignoring the triviality in d = 4 where the tree level marginal coupling constant 54 is 
actually irrelevant due to the radiative corrections. 
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renormalization group flow shown for d = 4 and 3 in Figs. 2 and 3, respec- 
tively. In those Figures the dimensionless coupling constants are displayed 
as the functions of the cutoff which is measured in the units of the initial 
cutoff value, k — > kA. 

In the case of theories with non-Gaussian asymptotically free couplings 
(the case of our model for d < 4), an excessive growth of these couplings 
in the UV regime may produce a non-perturbative situation in the infrared. 
We want to study this by comparing the values of the couplings at k = 0. To 
do so, we will adopt the convention that a model with positive renormalized 
mass square is non-perturbative in the vertex g n (j) n if the radiative correction 
0(g n ) to the self energy is stronger than the mass term, i.e. 

Uk) » 1 (41) 



((n-2)/2)!2(™- 2 )/ 2 m 2 (fc) 

where m 2 (k) = cj2{k)- In case this inequality were satisfied in the infrared 
(k = 0), this would mean a non-perturbative situation and the invalidity of 
renormalized perturbation expansion. 

We remark that we are asking here about the validity of the perturbative 
condition at the true, k = 0, infrared fixed point, where the RG flow ends. 
We are not considering for example the situation at the Wilson-Fisher IR 
fixed point in d < 4, which, as was explained in Section 4, is the IR limit of 
the UV scaling behavior (a fixed point which describes the behavior of the 
system at the end of the UV regime), and not the real IR fixed point if we 
ask for the behavior of the system at energy scales much lower than the mass 
m(k) (for example, in dimension two, we will see that at the crossover, or 
the end of the UV regime, the high order couplings start to take large values, 
while they are irrelevant couplings for the Wilson-Fisher IR fixed point; this 
is because the true RG flow separates from the extrapolation of the UV flow, 
as was explicitely shown in Fig. 1 in the case of d = 3). 

We now turn to a detailed analysis of the situation at dimensions four, 
three and two. Since the neglected higher order vertices may influence the 
evolution while we lower the cutoff, we have to address the problem of the 
system of coupled equations numerically. 

The evolution of a non-renormalizable coupling constant, g§ in four di- 
mensions, is shown in Fig. 2. The irrelevance is expressed by the indepen- 
dence of ge(k) on the initial value for k -C 10. The value given in the 

3 The different initial value for the non-renormalizable coupling constants may induce 
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leading order of the perturbation expansion,^] 0(gl), is reached at k ~ 0.3 
for not too large values of <7e(l)- Since the model is in the weak coupling 
regime the evolution is rather slow after arriving at this universal value if the 
cutoff is high enough to provide a long scaling regime. In our case the scale 
window 0.3 < k < 1 was insufficient and the plateau is reduced into a peak at 
k w 0.3 before the crossover. But the bringing of ge(l) close to the universal 
value creates a plateau even with this limited range of the scales as it can 
be seen for ge(l) = 9 • 10~ 8 . Such a scale independence is the ultimate goal 



of the improved action program |L2|]. As the initial value increases beyond 
the plateau level the coupling constant decreases in a monotonous manner 
as the cutoff is lowered. In order to check the critical exponent coming from 
the linearized blocking we need a g^(k) larger than the universal value since 
the latter originates in the nonlinear level, 0(gf). The evolution follows the 
linear relation g§ « k 2 for k < 1 indicated by the dashed line in the last plot 
of Fig. 2. The infrared scaling regime is < k < 0.1 where g§{k) tends to 
zero with k along the universal trajectory. The theory remains perturbative 
in d — 4 since it has no relevant non-Gaussian coupling constant. 

The three-dimensional renormalization group flow is depicted in Fig. 3. 
The asymptotic infrared scaling laws are rather simple, the super-renormal- 
izable coupling constants diverge, the renormalizable one n = 6 converges 
and the non-renormalizable ones tend to zero in the infrared, k — > 0. The 
ultraviolet scaling law, above k « 0.1 indicates the weak, radiative correction 
generated relevance of g^ for the given initial conditions, <? n (l)- The insensi- 
tivity on the initial condition cfe(l) for k < 0.3 seen in the last plot supports 
the irrelevance of g$. In fact, the evolution of gs{k) follows the linearized 
scaling law as long as its value is far form the universal 0(gf) value. It is 
worthwhile noting that the non-renormalizable coupling constants g n (k) al- 
ways develop a peak of sign (— l) 1+n / 2 around the crossover, k ~ 0.1. The 
appearance of the peak can be understood as the result of the increase of 
\g n (k) \ from zero as the cutoff is lowered in the ultraviolet scaling regime and 
the decrease in the infrared side of the crossover. 

Non-perturbative phenomena may arise at the low energy edge of the 
UV scaling regime due to the increase of the asymptotically free coupling 
constants, g± and g§ if the scaling regime is long enough and the initial value of 



a different overall scale factor. This effect is very weak in our case due to the smallness of 
the renormalizable coupling constants. 

4 Which is not applicable for the strongly coupled case ([38]) or (ff(J). 
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the coupling constants #4(1) and ge(l) are large enough. There is however no 
problem in finding a perturbative, asymptotically free theory in the infrared. 
The parameter with the highest energy dimension in the Lagrangian is #2 = 
m 2 for d > 2. Thus c/2 is the largest among the dimensionless coupling 
constants according to Eq. (p3|), and it dominates the action and renders 
the theory perturbative in the IR limit. The comparison of this conclusion 
with Eq. (^) reveals the necessity of treating the IR scaling laws properly 
in establishing the validity of the renormalized perturbation expansion. 

As d approaches 2 more and more coupling constants become super- 
renormalizable. The fastest increasing dimensionless non- Gaussian coupling 
constant during the decrease of the cutoff is g 4 . The critical exponent, the 
measure of the speed of the increase become degenerate for infinitely many 
coupling constants when d = 2. The specialty of the lower critical dimension 
is the existence of infinitely many super-renormalizable coupling constants, 
g n , with equal critical dimension. This degeneracy of the dimensions evolves 
the non-Gaussian pieces of the action with the same rate as the mass term 
on the tree level and the theories are not obviously perturbative any more. 
In other words, it remains for the radiative corrections, the last term in the 
right hand side of Eq. (^) to determine if the theory runs into weak- or 
strong-coupling regime at the infrared. 

The renormalization group equations were integrated out numerically in 
two dimensions with the initial conditions {g\ = —0.001, g\ = 0.01, 5^ = 
0, n > 4) to find the evolution of the coupling constants. The result, depicted 
in Fig. 4, shows a marked increase at the low energy end of the UV scaling 
regime. This increase originates from the asymptotically free evolution. The 
relevant behavior of a coupling constant is defined on the linearized level of 
the blocking, i.e. in the leading order of the perturbation expansion. This 
one-loop result can be obtained by replacing the running coupling constants 
in the beta functions by their initial values at the cutoff. The local potential 
obtained in the one-loop approximation is 

£4(0) = VM) + \( ^ log[ P 2 + TO)] = 

2 Jk<p<A ( 2 7TJ 

^{[A 2 + K'(0)]log[A 2 + TO)]- 

[k 2 + TO)] log[£; 2 + TO)] - A 2 + k 2 } , (42) 
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with 



W) = ™ 2 (A) + ^P0 2 - (43) 

The comparison of the numerical solution with the one-loop evolution is 
shown in Fig. 5. The one- loop formula ([42]) cannot be extended down to 
k = 0, because there will be a value of k such that k 2 +m 2 (A) = 0, since m 2 (A) 
is negative. But we only want to compare the results in the perturbative 
regime. So we have stopped the evolution in Fig. 5 at k ~ 0.3. One can also 
see an increase in the one-loop solution at small k which accumulates and 
drives the system non-perturbative at lower values of k. 

The numerical results of Fig. 4 show that the initial conditions (g\ = 
— 0.001, g l 4 = 0.01, = 0,n > 4) correspond to a non-perturbative system. 
Can we find initial conditions which yield perturbative dynamics? In order 
to answer this question the left hand side of the inequality (f4"T|) is plotted 
against the initial value for Q2 on Fig. 6, for the different couplings up to 
N = 20, at a value of g\ = 0.001. The result does not change qualitatively 
for different values of bare g^. It supports the general trend of having the 
systems more perturbative when the Gaussian part of the action is increased. 
The higher order coupling constants tend to grow faster but it seems that 
g n can be brought into the perturbative regime for sufficiently large initial 
mass square. At g\ = 0.01, for example, all the couplings are perturbative, 
and this perturbative character is more pronounced for the high couplings. 
However, the separation between the values of the l.h.s. ratio of Eq. ( f4~l"| ) at 
g\ = 0.01 is smaller as we go to higher couplings, and one can ask whether 
it has got a limiting value or the trend can be reversed for a sufficiently 
high order coupling. Indeed, Fig. 7 reveals that this happens for the n = 24 
coupling (N = 26), which suggests that one will have a non-perturbative 
situation also at this value of bare #2 going to a sufficiently high order coupling 
constant. The situation is the same, even stronger, below two dimensions: 
the existence of an infinite number of relevant couplings makes that one 
cannot assure perturbativity by looking to a finite number of couplings, no 
matter what the initial conditions are. We take this and similar other failures 
in finding a perturbative theory observed at different initial conditions as a 
strong numerical indication of the non-perturbative nature of any two and 
lower dimensional scalar field theory with polynomial interaction. 
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6 Lattice Wegner-Houghton equation 



The previous Sections dealt with the renormalization group flow at finite 
scales. We address now a different, asymptotic problem, the manner the sen- 
sitivity on the initial values of the irrelevant coupling constants is suppressed 
during the renormalization. This question is usually rendered trivial by the 
universality argument. But there are two reasons to suspect that such a rea- 
soning which is based on the linearization of the blocking relation might be 
oversimplified; both are related to an infinite set of operators. 

The reason motivating a more careful check of the universality, mentioned 
in the Introduction, is that the models at or below the lower critical dimension 
contain infinitely many relevant operators. It is not obvious whether the 
sum over the interaction vertices is always convergent enough to make the 
linearization of the blocking relation a reliable approximation. 

Another potential problem shows up if one changes infinitely many ir- 
relevant terms in the action by choosing another regulator. Let us compare 
the momentum space cutoff in the continuum with the lattice regularization. 
The propagator is a monotonic function of the momentum in the continuum. 
This is not the case on the lattice. In fact, the fermion doubling problem 
on the lattice [I3| results from the periodicity of the propagator in the first 
Brioulline zone, the appearance of 2 d — 1 new maxima in the propagators in 
the UV, non-universal regime. The existence of a maximum of the propagator 
in the UV regime contradicts an assumption of the studies of the continuum 
models, namely that the propagator decreases monotonically as j) — > oo, 
and renders the perturbation expansion non-universal for lattice fermionic 
models [14[]. There is no species doubling for bosons but their propagator 
remains periodic on the lattice and we find 2 d — 1 lattice extrema in the UV 
regime. The existence of these extrema is an evidence of the slowing down in 
the decrease of the propagator as the momentum approaches the boundary 
of the first Brioullin zone. This in turn indicates the weaker suppression of 
the high energy modes compared to the continuum regularization. Does this 
mean that the UV scaling laws are different in the continuum than on the 
lattice? We shall find an affirmative answer to this question but this result 
does not contradict the universality. 

Some words of caution are in order at this point. One would object the 
interpretation of the modification of the cutoff as the introduction of new 
irrelevant coupling constants by recalling that the theory ceases to be renor- 
malizable in the presence of the non-renormalizable (irrelevant) couplings. 
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The resolution of the apparent paradox is based on the difference between 
the ways the renormalization group is used in Statistical and High Energy 
Physics. We are interested in the dynamics close to the cutoff in Statisti- 
cal Physics and this is respected by the employment of the blocking which 
keeps the complete dynamics unchanged below the actual cutoff. The price 
of this precision is the appearance of the infinitely many irrelevant coupling 
constants in the action. We seek the dynamics at finite, fixed scales in High 
Energy Physics. Since the cutoff is sent to the infinity this boils down the 
problem of keeping the physics cutoff independent far from the cutoff only. 
The obvious gain of such an ease of the conditions is the freedom from the 
adjustment of the non-renormalizable parameters. Thus one can remove 
the cutoff when the non-renormalizable parameters are present in the action 
without any problem^ as long as the renormalization conditions are imposed 
far from the cutoff. 

The lattice regularization of the scalar model can be described by using 
the momentum space as the introduction of the non-renormalizable higher 
order derivative terms, 



{d^f — > -(sin^0) 2 (44) 




The cutoff dependence of the non-renormalizable coupling constants follows a 
tree-level relation arising from the Taylor expansion of the sine function. This 
is sufficient to establish convergent physics at finite scales when a — > |[14 



a claim to be verified in this Section numerically by means of the implemen- 
tation of the Wegner-Houghton scheme on the lattice. But this convergence 
can not rule out a modification of the scaling laws in the asymptotical UV 
regime. In fact, we shall find a new scaling regime between the region where 
the usual universal UV scaling is observed and the UV fixed point. The only 
effect the different adjustments of the non-renormalized coupling constants 
may leave on the finite scale physics can be comprised in an overall scale 
factor. 

It is rather straightforward to repeat the steps leading to Eq. (|TT|) on the 
lattice. It is shown in the Appendix that the only change required is the 



'Ignoring again the possibility of the triviality, the appearance of an UV Landau pole. 
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modification of the "solid angle" factor, Q(k): the lattice evolution equation 
(|7l|) is obtained from Eq. ([□]) by replacing Eq. (|1^) by Eq. (^). One recovers 
the continuum solid angle for d — 2, Q2 — 27r, in Eq. fl68|) as k — > 0, thus 
the WH equations agree in the IR limit. In fact, one sees numerically that 
the behavior of the evolution of the different coupling constants in the lattice 
RG is qualitatively the same as in the continuum case. But the question 
we are interested in is the relation between the regularizations in the UV, 
where the coupling constants are introduced, when the physics is the same 
at finite scales. It is shown in the Appendix that there is a natural relation 
between the cutoff's, A 2 = 8/a 2 , which matches the finite scale physics. We 
shall follow the renormalization group flow in terms of the coupling constants 
whose dimension is removed by the initial value of the cutoff, 

9n — > WA 2 , (45) 

in order to avoid the singularities at k = 0. 

Let us consider A0 4 lattice theory which can be studied either numerically 
or analytically and whose properties can be matched to those of the contin- 
uum theory by the adjustment of gi and g^. But the situation is more involved 
in two dimensions. The reason is again that there are infinitely many renor- 
malizable coupling constants and one cannot match the finite scale physics 
by adjusting g 2 and g^. This is demonstrated in Fig. 8 where the lattice 
model with the initial conditions m 2 (8) = #2(8) = —0.001, #4(8) = 0.01, and 
g n (8) = 0, n > 4 (where we have already made the rescaling Eq. ([45])) was 
evolved in the infrared direction. As the system reached k = k en< i the contin- 
uum WH equation was used to increase the cutoff. The result is a "perfect 
matching" of the models in the UV which gives the same low energy physics 
in the IR. As we can see in Fig. 8, the lattice A0 4 model in two dimensions 
is not the continuum A0 4 theory. It contains contributions of infinitely many 
renormalizable other coupling constants. Of course, numerically we had to 
truncate the equations at a certain coupling (here, at (9($ 22 )) but we checked 
that these "truncation effect" hardly influences the values of the low order 
coupling constants. 

We have taken for the parameter k 2 nd the value k 2 = 0.3 in Fig. 8, while 
the crossover is at k 2 T ~ 0.01. We had to use k 2 nd > k 2 T , because the high 
order couplings have very large values at the crossover which requires very 
fine discretization in the numerical resolution of our differential equations 
to ensure that the way back to the UV is done accurately. However, k 2 nd 
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should also be sufficiently small that the flow be universal there, in other 
words to make sure that the irrelevant lattice contributions are suppressed 
for k 2 < k 2 r . 

The choice k 2 nd > introduces an uncertainty in the matching. To assess 
it we repeated the "go-return" evolution described above and checked the 
discretization errors for k 2 nd = 0.8, 0.5 and 0.3. After then we took the 
appropriate bare parameters at the UV end points in both regularizations 
and followed the evolutions down to k — 0. The relative difference, 

. _ k CONT (0)- g ,r T (Q)l m 

- |<A ATT (0)| { ' 

is shown in Table 1 for the different coupling constants at k — 0. The 
smallness of the deviation assures that the IR behavior has practically been 
obtained with k 2 nd = 0.3, and one can trust the conclusions, the approach of 
the flow to a universal curve, extracted from Fig. 8. 

The first two plots in Fig. 8 show that the mass and the quartic coupling 
constant run parallel in the SUV region of the lattice regularization (see 
Appendix) and in the continuum. There is no convergence between the 
two regularizations in this unusual scaling regime, anticipated above. The 
approach to the universal curve starts for k 2 < 4, below the SUV regime only. 
The fact that the renormalization group flow converges to the universal one 
in the 2~ d -th part of the Brioullin zone only sets an unexpected high lower 
limit on the lattice size when the continuum limit is sought in numerical 
simulations. The higher order vertices seem to converge to the universal curve 
from the very beginning but the difference between the two regularizations is 
surprisingly large. The universal trajectory of the 4 model is reached later 
by the higher order vertices. This effect appears to be a counterpart of the 
non-perturbative features seen in Figs. 6, 7 and introduces a large uncertainty 
in identifying two-dimensional models in different regularizations. 



7 Conclusions 

The renormalization group flow of scalar models with polynomial interaction 
is considered in the first part of this paper by solving the Wegner-Houghton 
equation numerically in the local potential approximation for d — 2, 3 and 
4. 
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The numerical results showed in this paper suggest that the length of 
the UV scaling regime which is needed to generate non-perturbative dynam- 
ics in the infrared shrinks to zero as the number of the asymptotically free 
coupling constants tends to infinity. In other words, the A-parameter tends 
to the cutoff as the lower critical dimension is approached, d —>■ 2. Such a 
behavior limits considerably the values of the coupling constants for a pertur- 
bative system in dimension 3 and renders all two- and lower dimensional field 
theories with polynomial couplings non-perturbative. This makes the under- 
standing of the noncritical low dimensional condensed matter systems more 
involved. The one- dimensional models belong to first quantized quantum 
mechanics and our result is a manifestation of the failure of the convergence 
of the perturbation expansion for an anharmonic oscillator. 

Such a conclusion does not invalidate the well known results for two- 
dimensional systems, such as the applicability of the Bethe ansatz, bosoniza- 
tion and the availability of certain exact information for models with con- 
formal invariance. Instead, it makes the asymptotic state structure and the 
relation between the the dressed particles and the states created by the ap- 
plication of the field operator from the vacuum more involved. 

We found an interesting analogy between the infrared Landau pole of 
the confining four- dimensional Yang-Mills theories and the low dimensional 
scalar models which opens the possibility of an unexpected, nontrivial struc- 
ture in the asymptotic states in the low dimensional scalar models. Viewed 
with interest in particle physics our conclusion suggests that one can avoid 
the IR Landau pole by following the evolution of the non-renormalizable 
operators. 

How to find the non-renormalizable operators whose presence stabilises the 
theories at low energies? It is well known that massive Lagrangians generate 
trivial infrared scaling laws, i.e. the Gaussian mass term is the only relevant 
operator in the infrared scaling regime. This is because the fluctuations are 
exponentially suppressed beyond the correlation length so the evolution of 
the coupling constants slows down at the infrared side of the crossover. The 
theories with dimensional transmutation, i.e. dynamically generated scale 
parameter or infrared instability only can support non-perturbative dynamics 
in the IR scaling regime. Thus the operators sought should be relevant 
in the IR regime, their growth being fed by IR or collinear divergences. 
There are few known cases only where the low modes are controlled by non- 
renormalizable operators. These include the four fermion contact term in 
solids inducing the BCS transition [15], the higher order derivative terms in 
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the action which generate inhomogeneous vacuua [|TB|, the common element 
being the onset of a Bose-Einstein condensation ||. 

In the second part of the paper the infinitesimal renormalization group 
scheme is generalized for lattice regularization. The matching of the con- 
tinuum and the lattice regularizations is carried out numerically and the 
approach of the universal renormalization group flow is demonstrated for the 
two-dimensional </> 4 lattice model. This result suggests that the naive argu- 
ment for the universality, which is based on the linearization of the blocking 
relations remains valid in the presence of infinitely many relevant operators. 
Other potential troublemakers, the infinitely many higher order derivatives 
contained in the lattice kinetic energy do generate a new, "super UV" scal- 
ing regime but universality is restored at the IR end of the usual UV scaling 
regime. Another use of the lattice regulated version of the Wegner-Haughton 
equation is the estimate of the finite size effects in a non-perturbative man- 
ner. This provides a useful check of the thermodynamic limit of the numerical 
results obtained in general on small lattices. 



Appendix 



The details of the derivation of the Wegner-Houghton equation in the lattice 
regularization are given in this Appendix. Let us consider the scalar field 
theory regularized on a lattice of lattice spacing a = 1. We want to derive 
a WH equation similar to Eq. ([11]). We integrate over spherical shells in the 
momentum space for the continuum regularization, because the propagator 
has spherical symmetry. This is no longer the case on the lattice, where we 
have 



□ 



J^pJ, j3 M = 2sin 



TV 
2 ' 



(47) 



Let us see the surfaces of equal value of the lattice propagator in two dimen- 
sions by performing the following change of variables: 



(Px,Py) ► (p,0) 

4 sin 2 & + 4 sin 2 = p 2 
tan# = 

sm(p x /2) 



(48) 
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We can see in Fig. 9 the form of the curves of constant propagator for several 
values of p 2 . p can be identified as the "momentum scale" that runs from 
the cutoff at p 2 = A 2 = 8 to the IR p 2 = 0. It is also clear in that figure that 
the value p 2 = 4 separates two regimes, still in the ultraviolet region, that 
we could call super- UV (SUV), for 8 > p 2 > 4, and normal- UV regimes. For 
p 2 ~ 0, the lines are spheres, and our change of variables (f|^) reduces to the 
usual relation between cartesian and polar coordinates. 

The absolute value of the Jacobian of the transformation ( fj8|) is found to 

be 

J = Jp (49) 



(l - £ cos 2 fl) (l - £ sin 2 



where J p is the usual Jacobian for the polar change of variables, J p = p. 
The transformation fl4"8|) can be easily generalized to three and four dimen- 
sions; however, we can only treat analytically the integral that appears in 
the derivation of the WH equation in the d = 2 case. 

To derive the equivalent of the Wegner-Houghton equation (|TT|) i n the 
bidimensional lattice regularization we will start from Eq. (ffi), and calculate 
the trace by integrating in momentum space over a shell k — 5k < p < k, 
where p is the parameter we have introduced in Eq. (53), 

itrlogP + [/;'] = I|i!P J l 0g[ 4 s in 2 |+4 s m 2 | + C/i'] = 

^y 2 klog[k 2 + U^}m, (50) 

with 

n(jb) = / d0 — 8 (51) 

J V64- 16P + fc 4 sin 2 2fl 
We have to distinguish two different regimes in making the integration: (i) 
k 2 < 4. In this region the range of values for 9 is (0, 2ir). 

p2tt o 4 

n(k)= d9 =Vq(fc), (52) 

Jo V64 - IQk 2 + k 4 sin 2 26 ~t 
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where we have split the interval (0,2%) into four intervals (0,%/2), (%/2,%), 
etc. Let us consider Qi(k). With the change of variable x = tan 6 and the 
notation 



k 2 = 4 - k 2 
this integral can be brought into the form 

fii(jfe) 



2dx 



k Jo a/(x 2 + b 2 )(x 2 + b- 2 )' 



(53) 



(54) 



where b 2 = 4/k 2 . The integral in Eq. fl5"4"| ) is related to the elliptic integral of 
the first kind |TJ F((j),t), 



fti(Jfe) = F 



% k , -77 

-,-V8-k 2 
2' 4 



F(<p,t) 



da 



1 — t 2 sin a 



(55) 



The result is the same for the other integrals fli, i = 2,3, 4. But F(%/2, t) is 
a complete elliptic integral, which can be expressed in terms of the hyperge- 
ometric function JX?], [18 



F(a,(3; 7; z) 



r( 7 ) ^r(« + n)r(^ + n)z" 



r(«)rC9)^ T( 7 + n 



as 



0(^ = 2^^(1,1;!; 



A; 2 )fc 2 
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(56) 



(57) 



(ii) 4 < k 2 < 8. This is the SUV region. We split again the integral into 
four integrations in the corresponding quadrants. By using the same change 
of variables as above, we have to calculate 

f 4dr 
Qi(k) = / ; (58) 
J vM^ + 16) -4F(l + x 4 ) 

where 

k 2 = k 2 - 4. (59) 
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Special care is needed at the limits of integration (recall Fig. 9). It can be 
seen that the four integrals can be put together in the form 



2/k 



Adx 



k/2 ^/x 2 (k A + 16) - 4fc 2 (l + x 4 ) 
2 '~ k dx 



kJk/2 yj{x - k/2){2/k -x){x + k/2){x + 2/k) ' 



(60) 



This integral is related again [17] to an elliptic integral and an hypergeo metric 
function: 

(8 - k 2 ) 2 ' 



m = f 2 F 



'- --1 


= 16 V 1 




1 










_2'P 


k 2 ' 


V2' 


2 



(61) 



We would like to have a common expression for Q(k) for both cases (i) 
and (ii). From Eqs. ( p7|) and (|6l|) , we find in fact that the expressions differ 
in a factor 2 for k 2 = 4! The reason is that actually our integral is divergent 
at this point. From Eq. ([51]) we see that the divergent integral is 



n(k 2 = 4) 



2tt 



d# 



sin 201 



(62) 



(in fact, the hypergeometric function (56) converges in general only in the 
unit circle |^| < 1 ||17|| ). We will see, however, that this divergence is inte- 
grable during the RG evolution from k 2 = 8 to k 2 = 0, and therefore it has 
no physical significance. In order to have a consistent, single expression for 
the cases 8 < k 2 < 4 and 4 < k 2 < 0, we will make use of the following 
property of the hypergeometric functions [ I7 |: 



1 1 

2' 2 



l + z 



11 Az 
2'2' ' (l + z 2 ) 2 



< z < 1. 



(63) 



Let us consider the expression 
property (|63|) one finds, 



which is valid for k 2 < 4. Using the 



where we have set 



11 Az 
2 ' 2'2 ; ' (l + z 2 ) 2 



Az 



■2-(l + z)F[ 



(8 - k 2 )k 2 



2\2 



16 



(64) 



(65) 
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This equation has two solutions for X clS cL function of k: 



k '' S - k \ (66) 



8-k 2 k 

but only the first one is admissible for Eq. (0), because it gives z < 1 for 
k 2 < 4, while the second solution gives a value greater than 1 in this region. 
Now, if we define k = 8 — k 2 , then we have z = (8 — k )/k , and Eq. (J64|) 
becomes 

IK*) = 2,^(1,1; !;*>). (67) 

Comparing this last expression with the result (|61|) for the case 8 < k 2 < 4, 
we can finally write 

nm 167T n /l 1 , 2 \ . , 8-P 

n(k) = -y— F I -, -; 1; 2 ) , with z 



fc 2 V 2 2' y " P 

k 2 = k 2 if A; 2 = 8 ... 4, 
k' 2 = 8-k 2 if P = 4 . . . 0. 



and < _ _ " ' ' (68) 



The hypergeometric function F ~; 1; z 2 ) can be computed directly 
from its definition (|56|). One can obtain a high precision in the evaluation 
of the series with a reasonable number of terms (say, around 50) when z is 
not very close to 1, say, for 0.7 > z > 0. For 1 > z > 0.7 we have used the 



following alternative formula |]18 



n , „ . r(a + /3) ^ T(a + n)r(/? + n) r 

*■(«./*« + ft*) = TfTb^g („!)* W + D 

-V»(a + n) - + n) - log(l - z)}{\ - z) n , 
(|arg(l-*)| <7T,|l--z| < 1), (69) 



where [18 



,,(,) = iJ^M, (70) 

which gives a better convergence for the function F (|, |; 1; z 2 ) near z = 1 
because it is a series in the variable (1 — z 2 ). 
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In conclusion, our generalization of the WH equation (|TTD for a lattice 
regularization in two dimensions is 

where k is the parameter p of Eq. ([48]), and Q(k) is given by Eq. fl68f) . 
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FIGURE CAPTIONS 



Fig. 1: Numerical RG evolution for the (a) dimensionless g 4 coupling and 
(b) dimensional g 4 coupling in the d = 3 g 4 4 scalar theory. Black 
points result from the integration of the complete beta function with 
(?2 = 0.052, while white points show just the UV regime (see text). In 
(a), the IR limit of the UV regime is given by the fixed point (36). 

Fig. 2: Renormalization group flow for the dimensionless coupling constant 
ge(k) in d — 4. The initial conditions are <ib(l) = #4(1) — 0.01, together 
with the values of g 6 (l) = 0.0, 4 • 10~ 8 , 9 • 10" 8 , 10~ 7 for the different 
lines in plot (a), and ^(1) = 10~ 3 for plot (b). In this last plot the 
points of the numerical renormalization group flow are joined by the 
linearized scaling g 6 ss k 2 , shown by a dashed line. 

Fig. 3: (a) Renormalization group flow for the dimensionless coupling con- 
stants g n (k), n = 2,4,6 and 8 in d — 3. The initial conditions are 
g 2 (l) = «jf 4 (i) = 0.01, g n (l) = for n > 6. (b) Evolution of g 8 (k) 
with the same initial conditions that before except for gs(l) = 6 ■ 10~ 5 , 
together with the linearized scaling law g%{k) fa k, indicated by the 
dashed line. 

Fig. 4: Renormalization group flow for the dimensional coupling constants 
in d — 2. The initial conditions are #2(1) = <74(1) = 0.01, ^(1) = for 
n > 6. 

Fig. 5: Evolution of the dimensional coupling constants obtained in the 
one-loop approximation (squares) and numerically (circles) for d = 2. 

Fig. 6: The left hand side of Eq. (41) as a function of the bare mass square 
for n = 4, • • • , 20 in d = 2. 

Fig. 7: The same as Fig. 6 but now including the n = 24 coupling. 

Fig. 8: Matching of lattice and continuum regularizations; we start with 
lattice field theory at the UV, go down to the IR (circles) and then we 
come back to the UV with the continuum regularization (squares); see 
text. 

Fig. 9: Lines of equal value of the lattice propagator. To the inside, the 
lines drawn correspond to p 2 = 7.9, 7, 6, 5, 4, 3, 2, 0.5, 0.1. 
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TABLE CAPTIONS 

Table 1: Relative differences between the continuum and lattice coupling 
constants at k — 0, Eq. (46), after having matched the couplings at the 
UV using the parameter k^ nd . 
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TABLE 1 





= 0.8 


A; c 2 nd = 0.5 


fcend = 0.3 




0.0166 


0.0096 


0.0052 


A<? 4 


0.0028 


0.0011 


0.0002 


A<? 6 


0.0188 


0.0116 


0.0069 


A<? 8 


0.0199 


0.0135 


0.0092 


Afifio 


0.1793 


0.0992 


0.0494 


A^12 


0.0680 


0.0424 


0.0259 


A#14 


0.0522 


0.0340 


0.0221 


Agfie 


0.9300 


0.5558 


0.3226 


A#is 


0.2116 


0.1261 


0.0731 


A^20 


0.0584 


0.0389 


0.0259 
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FIGURE 1 




FIGURE 2 
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FIGURE 3 
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FIGURE 4 
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FIGURE 5 
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FIGURE 9 




